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A  GLOBAL  IONOSPHERIC  MODEL 


1.0  INTRODUCTION 

The  purpose  of  this  paper  is  to  describe  a  global  ionospheric  model  used  by  NRL  to 
analyze  the  performance  of  high-frequency  radar  systems.  The  basic  ionospheric  model  was 
developed  jointly  with  the  Institute  for  Telecommunication  Sciences,  and  was  primarily  based 
on  numerical  maps  of  ionospheric  indices.  These  numerical  maps  were  derived  from  reports  by 
a  network  of  ionospheric  vertical  sounding  stations.  Because  the  network  was  rather  sparse  in 
the  polar  region,  the  numerical  maps  do  not  accurately  portray  the  polar  ionosphere.  The  Air 
Force  Cambridge  Research  Laboratory  has  a  relatively  large  polar  ionosphere  data  base  which  it 
has  used  to  derive  a  set  of  polar  corrections  to  the  numerical  maps  of  the  F'2-layer  critical  fre¬ 
quency.  These  polar  corrections  are  described  in  Section  3.4  and  have  been  incorporated  into 
the  ionospheric  model  used  by  NRL.  Section  5  describes  some  of  the  problems  to  which  the 
ionospheric  model  has  been  applied. 

Section  6  provides  an  atlas  of  plasma  frequency  contour  maps.  This  atlas  is  a  graphical 
representation  of  the  ionospheric  model  and  will  be  useful  in  providing  an  intuitive  feeling  for 
the  model.  Some  simple  graphical  raytracing  may  also  be  carried  out  with  these  contour  maps. 

2.0  BACKGROUND 

For  many  years  numerous  organizations,  both  governmental  and  private,  have  been 
employing  the  HF  spectrum  to  communicate  point-to-point  between  long-distance  stations.  It 
was  recognized  early  that  HF  communication  systems  were  subject  to  marked  variations  in  per¬ 
formance,  and  it  was  hypothesized  that  most  of  these  variations  were  directly  related  to  changes 
occurring  in  the  ionosphere.  Considerable  effort  was  made  in  the  United  States,  as  well  as  in 
other  countries,  to  develop  research  teams  for  investigating  ionospheric  parameters  and  deter¬ 
mining  their  effect  on  the  nature  of  radio  waves  and  the  associated  reliability  of  HF  circuits. 
The  investigators  soon  realized  that  effective  operation  of  long-distance  HF  systems  increased 
in  proportion  to  the  ability  to  predict  variations  in  the  ionosphere,  since  such  an  ability  permit¬ 
ted  the  selection  of  optimum  frequencies,  antenna  systems,  and  other  circuit  parameters  that 
would  capitalize  on  ionospheric  variations.  With  the  encouragement  provided  by  these  findings, 
it  was  decided  that  more  raw  ionospheric  data  were  necessary  in  order  to  develop  models  that 
could  be  used  to  correctly  anticipate  ionospheric  conditions  affecting  HF  propagation.  World¬ 
wide  vertical-incidence  ionosondes  were  established  which  now  measure  values  of  parameters 
such  as  foE,  FoF\,  f„E^.  f„Fi,  and  h'F.  Worldwide  noise  measurement  records  were  started 
and  steps  wer  taken  to  record  observed  variations  in  signal  amplitude  over  various  HF  paths. 
The  results  of  this  research  established  that  ionized  regions  ranging  from  approximately  80  to 
600  km  above  the  earth’s  surface  provide  the  medium  of  transmission  for  electromagnetic  ener¬ 
gy  in  the  HF  spectrum  (3  to  30  MHz).  Furthermore,  most  variations  in  HF  system  perfor¬ 
mance  are  directly  related  to  changes  in  these  ionized  regions,  which  in  turn  are  affected  in  a 
complex  manner  by  solar-activity,  seasonal,  and  diurnal  variations,  as  well  as  latitude  and  longi¬ 
tude. 
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The  Radio  Propagation  Unit  of  the  U.S.  Army  Signal  Corps  provided  a  great  deal  of  infor¬ 
mation  and  guidance  in  1945  on  the  phenomena  of  HF  propagation  by  issuing  Technical  Report 
No.  6(1].  By  1948  a  treatise  of  ionospheric  radio  propagation  (2]  was  published  by  the  Central 
P.adio  Propagation  Laboratory  of  the  National  Bureau  of  Standards.  This  document  outlined 
the  state  of  the  art  in  predicting  expected  maximum  usable  frequencies  (MUF),  depicted  practi¬ 
cal  problems  of  ionospheric  absorption,  covered  in  detail  acceptable  methods  for  determining 
the  MUF  for  any  path  at  any  time,  and  took  into  account  the  various  possible  modes  of  propa¬ 
gation  by  applying  principles  which  were  found  to  work  in  practice.  The  model  used  to  make 
the  MUF  predictions  em,-loyed  the  "two-control-point"  method  and  assumed  the  ionosphere  to 
be  concentric,  with  reflection  occurring  only  from  the  regular  if  and  Fi  layers. 

In  1950  Laitinen  and  Haydon  of  the  U.S.  Army  Signal  Radio  Propagation  Agency  furth¬ 
ered  the  science  of  predicting  HF  system  performance  by  developing  empirical  ionospheric 
absorption  equations  and  combining  them  with  the  theoretical  ground  loss,  free-space  loss,  and 
antenna  gain  factors.  Thus,  expected  field  .strengths  could  be  anticipated  for  radio  signals 
reflecting  from  the  £and  Fi  regions,  considering  the  effect  of  solar  activity,  seasonal,  and  diru- 
nal  extremes.  These  findings  were  published  in  Technical  Report  No.  9  (3). 

The  accumulative  techniques  and  methods  presented  in  the  cited  literature  and  in  a 
number  of  other  studies  were  then  combined  to  establish  effective  manual  methods  for  predict¬ 
ing  the  expected  performance  of  HF  communication  systems;  however,  these  methods  were 
laborious  and  time-consuming  even  when  only  estimates  for  the  MUF  and  optimum  transmis¬ 
sion  frequency  (FOT)  were  needed.  To  alleviate  this  problem,  electronic  computer  routines 
were  developed  by  such  organizations  as  Stanford  Research  Institute  (1957)  [4),  Radio  Cor¬ 
poration  of  America  (1961),  and  the  Central  Radio  Propagation  Laboratory  (1961),  all  of  which 
were  based  upon  the  established  manual  prediction  methods.  The  CRPL  program  (5)  was  the 
first  computerized  technique  that  incorporated  a  numerical  coefficient  representation  of  the 
ionospheric  characteristics  (6].  However,  only  the  expected  MUF  and  FOT  were  predicted, 

In  1962  NBS  Report  7619  (7j  was  issued.  This  report  outlined  a  computer  routine  that 
utilized  the  then  most  recent  improvements  in  the  theory  of  performance  predictions,  combin¬ 
ing  the  more  predictable  ionospheric  characteristics  with  circuit  parameters  to  calculate  expected 
HF  system  performance:  MUF-FOT,  sys.em  loss,  reliability,  and  so  forth. 

In  1966  ESS  A  Technical  Report  IER-ITSA-1  (8]  was  published  with  an  improved  electron 
density  model.  In  this  model  the  electron  density  profile  along  the  path  was  assumed  to  be  ade¬ 
quately  represented  by  two  parabolic  layers;  that  is,  the  E  and  F^  layers.  The  height  of  max¬ 
imum  ionization,  thickness,  and  electron  density  were  derived  for  locations  near  the  points  of 
actual  reflection  along  the  path  instead  of  at  control  points  2000  km  from  each  end  of  the  path. 

Work  beyond  ITSA-1  was  continued  in  two  separate  paths,  one  for  communication 
analysis  and  predictions,  reported  in  ITS-78  [9],  and  another  for  analysis  and  prediction  of  the 
performance  of  over-the-horizon  (OTH)  radar  systems,  reported  in  i>'RL  Memorandum  Reports 
2226  [10]  and  2500  (!!].  The  electron  density  model  described  in  this  report  is  a  descendant  of 
the  routines  developed  for  predicting  the  performance  of  OTH  radar  systems. 
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3.0  IONOSPHERIC  INDICES  MODEL 

3.1  Introduction 

Prediction  of  ionospheric  indices  is  used  extensively  to  estimate  the  performance  of  long¬ 
distance,  HF  radio  systems  and  is  useful  in  the  design  of  earth-space  communication  systems. 

The  ionosphere  exhibits  considerable  statistical  variability.  If  the  minule-to-minute  varia¬ 
tions  within  the  hour  and  the  day-to-day  variations  within  the  month  are  averaged,  the  remain¬ 
ing  temporal  variations,  i.e.,  diurnal,  seasonal,  and  solar-cycle,  become  quite  well  behaved. 
These  remaining  variations  characterize  what  is  normally  referred  to  as  the  quiet  ionosphere 
because  the  percentage  of  disturbed  days  in  a  month  is  usually  relatively  small. 

It  is  the  purpose  of  this  section  to  review  what  indices  are  used  to  describe  the  quiet  iono¬ 
sphere;  tlie  next  section  will  describe  their  use  in  deriving  a  complete  electron  density  profile. 


3.2  The  Lower  Ionosphere 

Measurements 

Information  on  electron  densities  in  the  lower  ionosphere  (SO  to  90  km)  is  very  inade¬ 
quate,  primarily  as  a  result  of  limited  observations.  The  technical  problems  of  observations  are 
formidable,  and  the  interpretation  of  measurements  extremely  difficult. 

Predictions 

No  D-region  indices  are  included  in  the  present  prediction  model.  The  effects  of  the  D 
region  are  accounted  for  in  the  electron  density  model  by  extrapolating  from  the  E  region  using 
two  exponential  tails;  this  is  described  in  Section  4.2. 


3.3  The  E  Region 
Measurements 

A  large  volume  of  vertical  incidence,  ionosonde  data  has  been  collected  over  about  three 
solar  cycles,  and  many  features  of  the  E  region  are  therefore  well  known.  The  minimum  vir¬ 
tual  height  of  the  E  region  and  the  variation  of  maximum  electron  density  within  this  region  as 
a  function  of  time  and  geographic  location  arc  readily  obtained  from  the  ionograms.  The 
phenomenology  of  sporadic  E  has  been  investigated,  but  a  number  of  problems  remain 
unresolved. 

The  ^-region  indices  which  have  been  systematically  scaled  from  the  vertical-incidence 
ionosonde  records  include 
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Index  Definition 

f„E  The  critical  frequency  of  the  ordinary  component  of  the  flayer;  i.e., 

that  frequency  at  which  the  signal  from  the  ionosode  just  penetrates  the  E 
layer. 

h'  Es  The  minimum  virtual  height  of  the  sporadic- flayer,  measured  at  the  point 

where  the  trace  becomes  horizontal. 

f„E^  The  maximum  frequency  of  the  ordinary  component  of  sporadic  f  (f,). 

h'E  The  minimum  virtual  height  of  the  flayer,  measured  at  the  point  where 
the  trace  becomes  horizontal. 

J'hE^  The  blanketing  frequency,  i.e.,  the  lowest  ordinary  wave  .‘'requency  at  which 
the  f,  layer  begins  to  become  transparent,  usually  determined  from  the 
minimum  frequency  at  which  ordinary  wave  reflections  of  the  first  order  are 
observed  from  a  higher  layer. 

Prediciwiis 

The  regular  f  layer  is  predicted  using  three  indices:  the  montnly  median  value  of  critical 
frequency,  the  height  of  maximum  ionization  of  the  layer  and  the  ratio  of  /i,„f  to  semi- 

thickness  (t,„£).  Using  a  numerical  mapping  method,  Teftin  has  nroduced  numerical 
coefficients  representing  f„E  for  computer  applications  on  a  worldwide  basis.  They  were 
mapped  in  terms  of  latitude,  longitude,  and  universal  time  [12].  The  numerical  coefficients 
were  derived  from  measurements  taken  during  1958  (high  solar  activity)  and  1964  (low  solar 
activity).  A  linear  interpolation  procedure  was  used  between  the  representative  data  for  the 
high  (sunspot  number  =»  150)  and  low  (sunspot  number  =»  10)  solar  activity  periods  to  obtain 
/„f  estimates  at  all  other  phases  of  the  solar  cycle. 

An  examination  of  monthly  median  /»'f  observations  indicates  negligible  seasonal  or  geo¬ 
graphic  variation  in  the  minimum  virtual  height  of  f-region  ionization.  A  typical  value  is  110 
km.  When  the  lower  region  is  included  as  an  exponential  tail  of  the  flayer,  an  //,„f  of  110  km 
and  a  value  of  5.5  for  the  ratio  (//„,6V>'»|6)  are  used. 

The  median,  upper,  and  lower  deciles  of  /.f,  urc  available  from  numerical  maps  [13]. 
However,  because  of  its  probabilistic  nature,  sporadic  f  is  not  included  in  the  contour  maps. 
The  staiistics  of  the  sporadic  f  are  used  when  the  ionospheric  model  is  used  for  virtual  path 
tracing  [10] 

3.4  The  F  Region 

Mcasurmenis 

The  vertical-incidence  ionosonde  network,  with  its  long  series  of  measurements  over 
much  of  the  world,  provides  the  current  basis  for  f-region  prediction.s.  The  following  indices 
have  been  systematically  scaled  from  the  vertical  ionosonde  records  [14]  although  few  stations 
report  all  of  them. 
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Index 


Definition 


fnFi  The  critical  frequency  of  the  ordinary  component  of  the  Fi  layer,  i.e., 
that  frequency  at  which  the  signal  from  the  ionosonde  just  penetrates 
the  /'2  layer. 

M(3000)f2  The  ratio  of  the  maximum  useable  frequency  foi  a  distance  of  3000  km 
for  the  F^  layer  to  the  critical  frequency  of  the  layer. 

f„F\  The  critical  frequency  of  the  ordinary  component  of  the  F,  layer, 
i.e.,  that  frequency  at  which  the  signal  from  the  ionosonde  jus* 
penetrates  the  Fy  layer. 


h'F  The  minimum  virtual  height  of  the  flayer,  i.e.,  the  minimum  virt.ial 
height  of  the  night  flayer  and  the  day  fi  layer.  It  is  measured  at 
the  point  where  the  f  traces  become  horizontal.  (In  earlier  years,  the 
minimum  virtual  height  of  the  night  flayer  was  often  combined  with 
that  of  the  day  F-^  layer,  the  combined  tabulation  being  designated  li'FF2. 
In  these  cases,  the  minimum  virtual  height  of  the  F\  layer,  h'Fy,  was 
tabulated  separately.) 


h'Fi  The  minimum  virtual  height  of  the  F^  layer,  measured  at  the  point 
where  the  Fi  trace  becomes  horizontal. 


h.Fi  The  virtual  height  of  the  F^  layer  corresponding  to  a  frequency  J\ 

where  /  0.834  /„f2.  This  is  based  on  the  assumption  of  a  parabolic 

ionization  distribution,  which  is  usually  considered  justified  as  an 
approximation  near  the  maximum  of  the  F^  layer . 


Predictions 

The  F2  layer  is  described  by  three  indices:  monthly  median  value.s  of  critical  frequency 
/„f2,  height  of  maximum  ionization  //,„f2,  and  ratio  of  /i,„f2  to  semithickness  y,„F2.  The 
monthly  median  values  of  f„F2  and  the  M(3000)f2  factors  are  available  at  numerical  map 
coefiicients  in  terms  of  modified  magnetic  dip  angle,  longitude,  and  universal  lime  !15].  The 
data  for  the  mapping  were  from  the  years  1954  thiough  1958.  The  solar  activity  dependence  is 
accounted  for  by  a  linear  least  squares  fit  between  the  high  and  low  Zurich  sunspot  nuir.bers. 
The  height  of  maximum  ionization  is  found  by  first  determining  li,,r 2,  the  virtual  height  of  the 
F2  layer  at  0.834  /„f2.  A  geometric  formula  lipFi  accurate  to  within  6%  was  described  by 
Shimazaki  (16!: 


hpFi  =  -  176  +  1 490/ Af  (3000)  f  2.  (3-1) 

The  height  of  maximum  ionization  is  then  found  by  removing  the  retardation  caused  by  lower 
region  ionization; 


/>mF2  “  ^2  -  Ret. 


(3-2) 
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The  formulas  for  retardation  depend  upon  the  assumed  electron  density  profile  of  the 
lower  layers.  For  example,  the  D-E  ionization  can  be  a  parabolic  layer  with  an  exponential 
tail,  the  E-F  valley  a  linear  profile,  and  the  F\  layer,  if  present,  a  linear  or  parabolic  ledge. 
The  ratio  's  given  by  coefficients  in  terms  of  sun’s  zenith  angle  and  geomagnetic 

latitude  [8]. 

The  existence  of  the  /'|  layer  is  given  in  terms  of  a  maximum  solar  zenith  angle;  i.e.,  the 
F\  layer  exists  only  when  the  solar  zenith  angle  is  less  than  lonosonde  data  were 

analyzed  to  produce  a  map  of  Z„„x  (17); 

2m, IX  =  «K,i  +  *x.i  +  ^a^_2  +  ‘-'OS  (X),  (3-3) 

where  the  a^'s  and  b^s  are  coefficients  from  the  map  for  a  particular  month,  R  is  the  Zurich 
sunspot  number,  and  X  is  the  modified  magnetic  dip  angle  suggested  by  Rawer  [18].  When 
solar  zenith  angle  is  less  than  Zmax>  is  men  determined  by 

f„F]  =  4  /)(  /^  +  (^2  +  b2R)  cos  X  +  cos^;<  (3-4) 

where  the  o's  and  b's  are  coefficients  from  a  map  for  the  particular  month  and  x  is  the  solar 
zenith  angle.  The  height  of  maximum  ionization  {h,„F\)  is  given  by 

//„,F,  «=  165  4-  0.6428  X.  (3-5) 

where  x  is  given  in  degrees.  The  ratio  {h,„F\/y,„F\)  is  assumed  to  be  4.  The  values  thus 
derived  for  /;,„F|  and  y,„Fi  are  only  tentative.  If  the  height  of  the  F2  layer  at  f„F\  is  lower  than 
h,„F[,  then  h,„F\  and  y,„F^  are  adjusted. 

In  the  fall  of  1975,  NRL  received  from  Terrance  Elkins  at  the  Air  Force  Cambridge 
Research  Laboratory  a  computer  deck  containing  .software  whicii  could  be  used  to  "correct"  the 
f„F2  coefficients  provided  by  the  ITS  median  model  [6).  The  polar  correction  software  received 
from  AFCRL  has  been  adapted  for  use  with  the  NRL  ionospheric  model.  This  software  pro¬ 
vides  corrections  to  f„F2  as  a  function  of  magnetic  index  (A",,),  day  of  the  year,  universal  time, 
and  geographic  location.  The  polar  corrections  to  the  ITS  median  model  arc  described  in  detail 
in  AFCRL  Technical  Report  TR-75-0549  (191.  They  are  summarized  below. 

The  first  correction  applied  is  a  Kp  correction  where  6,  the  correction  factor,  is  a  function 
of  Kp  only.  The  Kp  correction  is  global.  All  of  the  other  correction  factors  are  applied  only  in 
the  polar  region.  The  equations  for  the  Kp  correction  are 

/of2  (corrected)  «/<,/’'2(I  -  8),  (3-6) 

8  =  0.025J^-0.I,  (3-7) 

where  the  term  X^  is  determined  by  the  following: 

If  Kp  <  0.3,  then  Xs  I; 


0.3 

< 

< 

1.3, 

Xn 

1.3 

< 

< 

2.3, 

Xn 

2.3 

< 

^P 

< 

3.3, 

Xs 

3.3 

< 

^P 

< 

4.3, 

Xn 

4.3 

< 

^P 

< 

6.3, 

Xn 

6.3 

< 

Kp, 

Xn 
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Thus,  a  Kp  value  of  3,  which  represents  a  normal  magnetic  activity  level,  produces  no  correc¬ 
tion.  Higher  magnetic  activity  reduces  /„/^2-  A  maximum  reduction  of  7%  occur?  when  Kp  > 
6.3. 

The  auroral  oval  correction  to  f^Fi  is  described  in  terms  of  magnetic  activity  (/f,,), 
corrected  geomagnetic  lime,  and  corrected  geomagnetic  latitude.  Software  for  determining 
corrected  geomagnetic  coordinates  was  also  obtained  from  AFCRL.  The  corrected  geomagnetic 
coordinate  system  used  was  described  by  Gustafsson  t2.0j.  For  the  auroral  oval,  the  correction 
factor  a  is  defined  such  that 


where 


fpF^  (corrected)  =  f„F2  (1  +  a), 


Here,  e  is  the  Napierian  base  and  k  is  defined  by 


(3-10) 


In  Eq.  (3-10),  X  is  the  corrected  geomagnetic  latitude  in  degrees  and  0  is  the  cquatorward 
boundary  of  the  auroral  oval  in  degrees.  X\  is  a  function  of  the  magnetic  activity  such  that 
4  <  A'l  <  6: 


AT,  -  7  -  K.,. 


(3.11) 


The  auroral  oval  correction  is  applied  if  the  corrected  geomagnetic  latitude  of  the  location  (X)  is 
greater  than  or  equal  to  </>.  This  equatorward  boundary  of  the  oval  is  described  by 

0  =  71.9-2.5  a;,- T.  (3-12) 

where  0  must  be  within  the  limits  0/^  <  0  <  0Ar.  with 

0^  =  68.9  - /(p  -  T.  (3-13) 

0,w  =  70.9  -  -  T.  (3-14) 

T  =  5.1  cos  (15  (T;.  -  D).  (3-15) 

Here.  7’.  is  the  corrected  geomagnetic  time  in  decimal  hours  (angle  is  in  degrees). 

The  trough  correction  is  the  next  to  be  applied.  This  is  used  only  equatorward  of  the 

aurora!  oval  and  only  when  the  solar  zenith  angle  x  >s  greater  than  90°  (nighttime).  The  trough 

corrcciion  is  not  applied  for  corrected  geomagnetic  times  in  the  range  6  <  T,  <  18.  For  the 
trough  correction  the  form  of  the  correction  factor  a  (applied  as  in  Eq.  (3  8))  depends  on  the 
value  of  k  (defined  in  Eq.  (3- 10)): 


n  (for  k  >  1), 


o  =  i.6487  like  ^  (for  k  <  1), 


where 


r,  =  0.2(1  +  cos(27rZ)/365))c‘-T^/'^' 

and  D  is  the  day  of  the  year.  The  term  y  is  defined  by  Eqs.  (3-19)  and  (3-.i0): 


(3-16) 


(3-17) 


(3-18) 


'Ur''* . 
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y  =  7;  -3  (forO  <  r,  <  6),  (3-19) 

y  =  27  -  T,  (for  18  <  7;  <  24).  (3-20) 

Note  that  only  these  two  intervals  need  be  defined  since  the  trough  correction  is  not  applied  in 
the  interval  6  <  T,  <18. 


In  the  interval  90  <  x  (solar  zenith  angle)  <  94.6,  the  factor  i\  is  reduced  by 


/i  =  I, 


Y-90 

4.6 


to  provide  a  smooth  transition  between  the  normal  ionosphere  and  the  trough. 


(5-21) 


3.5  Nummary 

The  ionospheric  .ndicc£  used  in  this  model  are  a  minimum  selection  froni  those  available. 
The  li'F,  Fi  maps  are  not  used  because  the  frequency  at  which  they  appear  is  not  available,  so 
it  is  not  possible  to  properly  adjust  for  retardation.  In  fact,  if  the  retardation  is  not  correctly 
accounted  for,  the  bottom  of  the  F  layer  occasionally  will  bo  calculated  as  higher  than  hmFi. 
Therefore,  the  map  of  the  ratio  (h,„F2/y,„F2)  is  used  to  give  consistent  results.  The  maps  of 
f„F2  and  ratio  /t,„ Fj/j,,, F2  are  from  the  same  data  base  and  are  statistically  consistent.  Using 
the  ratio  guarantees  that  y,„F2  is  positive. 

4.0  ELECTRON  DENS»  rv'  PROFILE  MODEL 

4.1  Introduction 

Section  3  described  how  a  set  of  indices  describing  the  vertical  electron  density  profile  is 
generated.  This  section  will  de.scribe  the  procedure  used  to  generate  a  vertical  electron  density 
profile  froru  this  set  of  indices.  The  electron  density  in  terms  of  the  plasma  frequency  squared 
is  given  by 

A' -  1.24  X  lO'Vw^  (4-1) 

where  N  is  the  number  of  electrons  per  cubic  meter  and  ,/)v  is  the  plasma  frequency  in  MHz. 
Just  prior  to  exiting  from  the  routine,  the  profile  is  converted  back  to  plasma  frequency  vs 
height  for  use  in  the  contour  mapping  routines. 

4.2  D  Region 

The  profile  is  generated  in  two  steps.  First,  the  coefficients  which  describe  the  various 
segments  are  calculated  and  then  the  electron  density  profile  is  generated  from  these 
coeflicients.  The  lower  D  region  is  considered  first.  The  lower  D  layer  (40  to  65  km)  is  given 
by 

/n  ^  Oi)  =  /w  (40)  pO  (4-2) 


where  y)v(40)  is  the  plasma  frequency  at  40  km,  which  is  assumed  to  be  2.01  x  10“^  MHz,  and 
//  is  the  height.  This  corresponds  to  an  electron  density  of  5  electrons  per  cubic  centimeter. 
This  model  of  the  lower  D  region  was  suggested  by  Nestorov  (21). 
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The  expression  used  to  describe  the  upper  D  region  (65  to  98  km)  is  an  adaptation  of  the 
exponential  model  suggested  by  Nestorov: 

/a^(/»)  =/Ar^(65)  (4.3) 

Rather  than  tie  the  upper  exponential  to  the  solar  zenith  angle  as  Nestorov  did,  the  exponential 
coefficient  k  is  chosen  to  merge  the  exponential  upper  D  region  with  the  £-region  model.  Thus 
the  upper  Z>-regiDn  model  follows  the  diurnal  and  seasonal  variations  of  the  ZT-region  model. 
The  expression  for  k  is 

k  =  \n \fs Hn)lfN H65) )/(98  -  65) ,  (4-4) 

where  //vf(65!  is  the  plasma  frequency  at  65  km  from  the  lower  Z)-region  model.  /a/(65)  is  8.98 
X  10"^  MHz.  This  corresponds  to  about  100  electrons  per  cubic  centimeter.  7^1(98)  is  found 
by  evaluating  the  fi-region  parabola  at  98  km.  The  slopes  of  the  equations  used  for  the  D  layer 
are  not  continuous  at  the  two  merge  points  (65  and  98  km).  This  produces  cusps  at  these 
points  in  the  virtual  height  profile.  However,  the  electron  density  is  sufficiently  low  at  these 
points  so  that  the  cusps  do  not  significantly  affect  the  direction  of  rays  in  the  high-frequency 
range. 


4.3  E  Region 


The  ionospheric  indices  used  are  f„E,  median  (ordinary-ray)  critical  frequency  of  the  E 
layer,  h,„E,  height  of  maximum  ionization  of  the  flayer  (110  km);  and  >>„£',  the  f-layer 
semithickness  (20  km).  The  electron  density  (proportional  to  plasma  frequency  squared)  is 
modeled  as  a  parabola: 


/AfO 


li„.  E  ~  h 
ymE 


(4-5) 


where  is  the  plasma  frequency  at  a  height  li. 


4.4  E  -  Fa  Valley 


Normally  there  is  a  valley  in  the  distribution  of  electron  density  between  the  E  and 
regions.  Only  the  total  density  in  this  valley  is  modeled,  not  the  shape.  In  this  area  between 
the  E  and  Fi  regions  the  electron  density  is  modeled  by  a  straight  line  between  a  point  on  the 
top  side  of  the  f-region  parabola  and  a  point  on  the  lower  side  of  the  fa'fegion  parabola. 
These  two  points  are  both  defined  in  terms  of  the  f„E.  The  point  on  the  top  side  of  the  E- 
region  parabola  is  at  0.8516  f„E.  The  point  on  the  bottom  side  of  the  fi-rcgion  parabola  is  at 
0.98  f„E.  The  height  li„  of  the  upper  point  is 


I'u  =  llml'2  -  ymf^2 


0.98  f„E 
LF 


(4-6) 


The  height  /;/  of  the  lower  point  is 

/)/,  =  h,„E  +  y„,E 


/  0.8516  J„E 

y  ‘  ~  r..E 


(4-7) 


The  constants  0.8516  and  0.98  have  been  chosen  to  represent  as  nearly  as  possible  measured 
depths  of  the  valley. 
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4.5  F|  Region 


The  F\  layer  is  described  by  three  parameters:  the  critical  frequency  (y„/’i),  the  height  of 
maximum  ionization  (h, „/'{),  and  the  semithickriess  0',„/'t)-  layer  nay  be  either  linear 

or  parabolic.  If  linear, 

(4-8) 

where  5|  is  the  slope  defined  by 


The  term  is  defined  by 


ys 


T.V  111  -  limP\  +T,„F, 


or 


Ts-  =  I . 

whichever  is  larger.  The  term  hi  is  height  in  the  Fi  layer  at  f„F\. 


Ill  “  hmPl-ymFl  1/  I  - 


/ 

i  — 

1 

y;./‘? 

If  the  F\  layer  is  parabolic, 


Sn^Ui)  “  (y;./'i)^ 


\  />»,/=’!  -  ll 


\  ymP\ 


(4-9) 


(4-10) 


(4-11) 


(4-12) 


The  choice  o'’  a  linear  or  parabolic  shape  to  the  electron  density  of  the  F|  layer  is  made  by  com¬ 
paring  the  luight  of  maximum  ionization  of  the  f|  layer  ili,„F\)  to  the  /•2-layer  height  at  f„F^ 
(111  as  defined  by  Eq.  (4-1 1)).  If  hi  is  higher  than  h„,Ft,  then  the  parabolic  layer,  Eq.  (4-12),  is 
used  for  the  F[  layer,  if  hi  is  not  higher  than  //.„/’i,  the  slope  S]  defined  by  Eq.  (4-8)  is  com¬ 
pared  with  the  slope  Fi  of  the  Fi  layer  at  the  point  hi  (frequency  is  /„/')): 


52 


2if„Fi)^  ih„,Fi  -  hi) 
^yml^l)^ 


(4-13) 


if  the  dififercnce  (5i  -  Si)  is  positive,  the  linear  Fi  layer  is  used.  If  the  difference  is  negative, 
the  parabolic  layer  is  used. 


4.6  F2  Region 

The  ionospheric  parameters  used  for  this  region  are  the  median  critical  frequency  (J„Fi), 
the  height  of  maximum  ionization  Ui,„Fi),  and  the  semithickness  iy,„Fi),  The  nose  of  the  Fi 
layer  is  described  by 


fsHh)  =  (f„Fi)^ 


h„,Fi  -  h 

2 

(4-14) 
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The  top  side  of  the  region  is  modeled  by  merging  a  single  exponential  to  the  top  of  the  Fi 
layer  at  a  height  /;,  determined  by  the  relation 


+  0.25  y„,F2. 

Above  height  li,  the  plasma  frequency  /(jih)  is  described  by  the  equation 


/nHIi) 


k(f,F2)^e 


Iz'j. 

~  k 


where 

k  =  1.875  y„,F^ 

This  top-side  model  is  an  adaptation  of  one  developed  by  Bent,  et  al.  (22l. 


(4-15) 


(4-16) 


(4-17) 


The  value  of  k  is  chosen  to  make  the  Fj-region  parabola  (Eq.  (4-14)),  the  top-side 
exponential  vEq.  (4-16)),  and  their  derivatives  continuous  at  /;,,  the  merge  point.  This  is  done 
to  keep  from  causing  a  cusp  at  the  merge  point  in  the  virtual  height  profile. 


5.0  APPLICATIONS 

5.1  Anomalous  Propagation  Diagnostics 

Plasma  frequency  contour  maps  can  be  used  to  predict  when  anomalous  propagation  con¬ 
ditions  may  occur.  For  example,  the  contour  map  shown  in  Fig.  1  exhibits  a  definite  positive 
tilt  from  the  origin  out  to  about  2700  n.mi.  Line  of  sight  at  zero  elevation  angle  is  shown  on 
the  contour  map  with  a  broken  line.  Since  low-elevation-angle  rays  are  more  likely  to  become 
trapped,  lilts  in  the  immediate  vicinity  of  this  broken  line  are  critical  for  trapping  during  the 
first  refraction.  An  estimate  of  the  operating  frequency  which  will  produce  elevated  modes  for 
a  particular  contour  map  can  be  obtained  by  finding  the  highest  plasma  frequency  encountered 
by  the  zero  elevation-angle  line  of  sight  and  using  the  well  known  secant  law  to  estimate  the 
equivalent  oblique  frequency  /„: 


/«  =  /.•  sec  (f>. 


(5-1) 


where  /,,  is  the  vertical  plasma  frequency  and  i,^'  the  angle  between  the  ray  line  of  sight  and 
the  zenith  al  the  true  height  of  reflection.  When  the  elevation  angle  is  zero,  the  angle  </»  may 
be  found  by  solving 

Ji 


whc>. 
sin^  ({} 


earth  radius  and  li  is  the  true  height  of  reflection. 


cos  i) 


1,  Eq.  (5-1)  may  be  rewritten 

/,.  =  kf^, 

where  k  is  defined  by  the  expression 


(5-2) 

Making  use  of  the  identity 


(5-3) 
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(5-4) 


11 


PLR5Mfl  FREQUENCY  CONTOURS 

LDCFiL  TIME 


WM-iHSrSH 


Example  of  conlour  map  which  indicates  that  an  elevated  mode  may  be  foond  by  raytracing.  (Broken 
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Following  the  broken  line  in  Fig.  1  upward  from  the  origin,  we  find  that  the  vertical  plasma  fre¬ 
quency  increases  until  it  reaches  a  maximum  of  about  5.5  MHz  at  280  km.  Equation  (5-3)  may 
be  used  to  predict  the  critical  frequency  for  the  zero-elevation-angle  ray.  In  this  case  we  would 
estimate  that  at  operating  frequencies  above  19.2  MHz  all  rays  would  penetrate  and  that  at 
operating  frequencies  below  19.2  MHz  the  low-angle  rays  would  be  refracted.  Because  of  the 
positive  lilt  in  that  region  of  the  ionosphere  where  these  rays  are  refracted,  the  low-angle  rays 
will  probably  be  tilted  into  an  elevated  propagation  mode. 


The  trajectory  of  the  zero-elevation-angie  ray  may  be  predicted  in  greater  detail  by  using  a 
raytracing  program  with  a  numerical  representation  of  the  contour  maps.  The  contour  map 
shown  in  Fig.  1  was  used  with  the  raytracing  routine  developed  by  Jones  and  Stephenson  [231 
to  trace  the  low-elevation-angle  rays.  Interface  between  the  raytracing  program  and  the 
plasma-frequency  contour  map  was  accomplished  by  an  adaptation  of  an  interpolation  routine 
written  by  ARGON  Corporation  for  AFCRL  [24].  This  routine  supplies  plasma  frequency  and 
the  required  spacial  derivatives  when  a  grid  of  plasma-frequency  vertical  profiles  is  input.  For 
an  operating  frequency  of  19  MHz  and  the  ionosphere  shown  in  Fig.  1,  the  raytracing  program 
predicts  that  the  zero-elevation-angle  ray  will  reach  a  maximum  altitude  of  276  km  at  1170 
n.mi.  before  being  refracted  back  toward  the  earth.  During  the  first  refraction  the  ray  will  be 
tilted  sufficiently  to  miss  the  earth.  The  point  of  closest  approach  of  the  ray  to  the  earth 
between  the  first  and  second  refractions  is  113  km  at  2000  n.mi.  Figure  2  shows  graphically  the 
trajectory  of  rays  between  0°  and  10®  elevation. 


The  point  of  all  this  is  not  that  the  propagation  mode  predicted  will  be  duplicated  precisely 
but  that  elevated  modes  of  the  type  predicted  by  the  raytracing  routine  may  be  expected  for  the 
combination  of  operating  conditions  examined. 


When  used  for  raylracing,  the  ionospheric  model  must  be  used  with  a  magnetic  field 
model.  Typically,  the  earth-centered  dipole  model  provided  with  the  Jones-Slephenson  raylrac¬ 
ing  program  is  employed.  In  this  model  the  gyrofrequency  fn  is  given  by 


///  =■  /// 


\i 


R 

R  +  h 


(1+3  cos^O)®-^, 


(5-5) 


where  ///  is  the  gyrofrequency  at  the  equator  on  the  ground  (typically  0.8  MHz);  R  is  the  radius 
of  the  earth;  h  is  the  height  above  the  earth;  and  9  is  the  geomagnetic  colatitude.  The  magnetic 
dip  angle  I  is  given  by 

tan  /  =>  2  cot  0.  (5-6) 


Occasionally  the  calculation  of  absorption  in  conjunction  with  raytracing  is  desired.  The 
Jones-Stephenson  raytracing  routine  provides  this  option.  The  software  package  supplied  with 
the  basic  raytracing  routine  includes  several  collision  frequency  models.  Typically  the  simple 
exponential  profile  (EXPZ)  is  used  where  the  collision  frequency  v  is  defined  by 


V  = 


(5-7) 


where  /;  is  height  above  the  ground;  //q  is  the  reference  height  (typically  70  km);  i^o  'S  the  colli¬ 
sion  frequency  at  the  reference  height  (typically  8  x  10®  collisions  per  second);  and  a  is  the 
exponential  decay  coefficient  (typically  0.16). 
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Comparisons  have  been  made  between  the  loss  predicted  by  the  raytracing  routine  and 
that  predicted  by  virtual  path  tracing.  For  example,  for  a  test  case  where  the  sunspot  number 
was  28  and  the  solar  zenith  angle  was  63°,  the  virtual  path  tracing  model  (10]  predicted  39-dB 
one-way  nondeviative  (Z)  region)  absorption  loss  for  a  1°  elevation-angle  ray  at  5  MHz.  This 
same  program  also  calculates  deviative  absorption  which  occurs  as  the  ray  penetrates  deeper 
into  the  ionosphere  and  deviates  from  a  linc-of-sight  path.  For  the  test  case,  the  deviative 
absorption  was  6  dB,  making  a  total  absorption  loss  of  45  dB.  Using  the  ionospheric  model 
described  in  this  report  with  the  Jones-Stephenson  raytracing  routine,  the  total  absorption  loss 
calculated  for  the  test  case  was  41  dB.  This  level  of  agreement  between  the  two  programs  gives 
some  confidence  in  the  .Cl-region  models  used.  However,  there  remains  one  major  defect  in 
using  the  Jones-Stephenson  raytracing  routine  with  this  ionospheric  model  for  loss  calculations. 
Because  of  the  probabilistic  nature  of  the  sporadic  E,  this  layer  is  not  handled  correctly  by  the 
raytracing  routine,  which  is  basically  deterministic.  In  the  virtual  path  tracing  program  a 
separate  loss  calculation  is  made  for  the  effect  of  sporadic  E.  This  loss  is  termed  obscuration 
loss  and  is  applied  to  F-layer  modes  in  the  virtual  path  tracing  program.  The  addition  of  this 
feature  to  our  version  of  the  Jones-Stephenson  raytracing  routine  is  planned. 

5.2  High-Frcqucncy  Radar  System  Performance  Prediction 

The  global  ionospheric  model  may  be  used  to  predict  the  performance  of  high-frequency 
OTH  radar  systems.  This  requires  that  a  model  of  the  radar  be  coupled  to  the  ionospheric 
model.  NRL  Memorandum  Reports  2226  (10)  and  2500  [11)  describe  how  an  earlier  version  of 
the  ionospheric  model  was  used  for  this  purpose.  This  procedure  makes  use  of  virtual  path 
tracing  and  therefoi?.  will  not  describe  as  completely  the  propagation  modes  associated  with  any 
specific  set  of  conditio.ns.  This  is  especially  true  during  transition  when  anomalous  propagation 
modes  exist,  f'or  these  modes,  the  ionospheric  model  should  be  used  with  the  Jones- 
Stephenson  raytracing  program  as  described  in  the  previous  section. 

6,0  ATLAS  OF  PLASMA  FREQUENCY  CONTOUR  MAPS 

To  provide  a  graphic  picture  of  how  the  plasma  frequency  varies  as  a  function  of  location, 
season,  time  of  day,  and  magnetic  activity,  a  set  of  plasma  frequency  contour  maps  are  it  cluded 
(Figs.  3-38).  All  of  these  contour  maps  begin  at  the  equator,  follow  the  69'’W  meridian  north 
through  the  north  geographic  pole,  and  then  south  along  the  111°E  meridian  back  to  the  equa¬ 
tor.  Thus  each  contour  map  provides  a  cross  section  of  the  entire  northern  hemisphere.  Maps 
are  provided  for  three  seasons,  equinox'  (represented  by  March),  summer  (represented  by 
June),  and  winter  (represented  by  December).  For  each  season,  contours  were  mapped  for 
0500  UT,  1100  UT,  1700  UT,  and  2300  UT.  These  times  correspond  roughly  to  midnight, 
morning,  noon,  and  evening  local  time  along  the  69®  meridian.  F'or  each  combination  of  time 
and  season,  contours  were  drawn  without  the  polar  corrections,  and  with  the  polar  corrections 
for  two  levels  of  magnetic  activity  (A^  •=  3  and  =■  7). 


Fig.  4  —  Mar-  ■■  map.  0500  UT.  polar  corrections  (A' 


Fig.  5  —  March  map.  0.^00  UT.  polar  correclions  (K, 


PLRSMB  FREQUENCY  CONTOURS 

LDCRL  TIME 


NRL  REPORT  8321 


wM-ii^raH 


19 


Fig,  6  —  March  map.  1 100  UT.  no  polar  correciions 


March  map.  1 100  UT,  polar  correciions  (/T, 
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P'lg.  8  —  March  map.  1 100  UT.  polar  correciions  (A' 
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March  map,  1700  UT.  no  polar  correciions 
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Fig.  M  "■*  tVfarch  map,  7300  UT,  polar  correciions  iK 
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Fig.  20  —  June  map.  1100  UT.  polar  correclions' ( AT, 
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Fig.  21  —  June  map.  1700  UT.  no  polar  corrections 
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Fig.  23  —  June  msp.  1700  UT.  poKir  eorrcclions  (A', 
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Fig.  27  —  December  map,  0500  UT.  no  polar  correclions 
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Fig.  29  —  December  map,  0500  UT.  polar  corrections  ( K, 
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Kig.  30  —  December  map.  1 100  UT.  no  polar  correclions 
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December  map.  i700  UT.  no  polar  corrections 


PLfiSMfl  FREQUENCY  CONTOURS 

LDCOL  TIME 


Fig.  34  —  December  map.  POO  UT,  polar  corrections  (A' 
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Fig.  35  —  December  map.  1700  UT.  polar  corrections  iK, 
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Fig.,  38  —  December  map.  2300  UT,  polar  corrections  ( K, 
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